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1. Introduction 

The incompressible Navier-Stokes equation 

dtU+{uyu)-vW^u= --Wp (1.1) 

P 

provides the basis for the description of laminar flow. Here u = u{x,t) is the velocity, p = p{x,t) the 
pressure, and v the kinematic viscosity. In the limit where v ^ 0+ however, the flow is known to 
become highly irregular, turbulence ensues, and one may ask if ( |1.1[ ) still captures the full dynam- 
ics. For such a state it is certainly necessary to consider a statistical theory where u{x, t) is a random 
field in space and time and ( |1.1| ) is replaced with an infinite hierarchy of dynamical equations that 
relate the different orders of correlation functions. Usually one has to rely on certain closure as- 
sumptions to close this set of equations (see e.g. [[I]]). One may take another route however, where 
one starts from reasonable assumptions on the symmetries of the problem. In particular, assuming 
statistical homogeneity, isotropy, and scale-invariance, simple dimensional analysis yields a tight 
prediction for the moments of velocity differences (structure functions) 

(A,m)« oc , Ch=«/3, (1-2) 

where Ar = (^(.x + r) — and the bar ••• denotes spatial averaging. The universality con- 

jecture then states that this scaling behavior should hold far from the boundaries and indepen- 
dent of the mechanisms that generate the flow. Nevertheless, both experiment and direct numerical 



simulations of ( |1.1| ) indicate a violation of this scaling behavior for high-order moments (see e.g. 
[^, 0]). In terms of symmetries this corresponds to the breaking of scale-invariance. Highly 
erratic, intermittent structures give the dominant contribution to these moments, the universal sta- 
tistical properties of which are still largely unknown. What is the nature of these structures, and is 
it possible to understand their properties from first principles? 

Here, the random-force-driven Burgers equation (see [|6[ ^ ||] for a review) is taken as a one- 
dimensional model 

dtU + udxU — vd^u= f (1-3) 

of the Navier-Stokes equation. We artificially generate a turbulent state by driving the system by a 
self-similar forcing that is white-in-time 

{f{k,t)f{k',t')) =Do\k\f^5{k + k')5{t-t') . (1.4) 

The brackets (• • • ) denote ensemble averaging, Dq is a dimensionful parameter, and /3 measures 
the relative strength of the forcing at different scales. For large, negative values of j8 the forcing 
effectively acts at large scales ~ L. On the other hand kinematic viscosity v provides a dissipation 
scale T] , and for v — > 0+ these two characteristic scales separate. In particular, for j3 = — 1 the 
interplay of the stochastic forcing and advective term leads to a Kolmogorov energy spectrum 
the intermediate range of scales, reminiscent of Navier-Stokes turbulence [Q, p^ ]. 
The physical picture that one may associate with this scenario is the appearance of shocks with 
a finite disspative width (see Fig.|l]). These structures give the dominant contribution to the high 



order moments of velocity differences, and leads to a strong form of intermittency 1 1 1, 12] where 



(|AM|")ocr, n>3. (1.5) 
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Figure 1: Velocity profile u{x) from a simulation on a 254 x 1024 (space x time) lattice, where x is taken 
in units of the spatial lattice size L. 



In view of the well-established anomalous scaling behavior of Burgers turbulence [13] and the 
physical picture of the underlying mechanisms for intermittency [I^, |ll], |l^], Burgers equation 
(1.3) provides an ideal benchmark setting to test new analytical and numerical methods for Navier- 
Stokes turbulence. 



2. Functional Integral 

The functional integral representation for the random-force-driven Burgers equation is ob- 
tained via the Martin-Siggia-Rose formalism by means of an auxiliary response field /X [14, 15, 16, 
17]. We have the partition function 

Z = j d[u]d[lx\&x^{-S[u,ix\} , (2.1) 

where the action 5[m,/i] is given by 

S = —ijdtdxii{dtU + udxU — vd^u)-\-—jdtdxdy}x{x,t)D{x — y)}x{y,t) . (2.2) 



Here, D{x — y) is the spatial part of the two-point correlation function (|E4|). Notice, that in this 
form the action does not satisfy positivity. To obtain a Gibbs measure that can be sampled by a 
Markov chain Monte Carlo (MCMC) algorithm we integrate out the auxiliary field. This leaves us 
with the probability density functional 

P[u] = exp I ~ 2 y dtdxdy{dtU + ud^u — V d^u)D^^ {x — y){d,u + ud^u — v5^m)| , (2.3) 
which is the starting point for our investigations. 



3. Lattice Theory 

The theory is defined by placing the field u{x,t) on the sites of a regular space-time lattice 
A, i.e. {x,t) € A. This way, we impose a UV cutoff that eliminates the details of those processes 
occurring deep in the dissipative regime. Then, the measure is given by d[u] n(v,f)eA'^"(-^)0 
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and the action in (2.3) needs to be discretized appropriately. We replace the dynamics (|L^ with a 
finite-difference equation with backward-time discretization 



dfU + udjcU - {u(t) — u{t — e)) + u{t — e)djcu{t — e) 



(3.1) 



where e is the lattice spacing in time direction. This ensures the correct dynamics in the continuum 



limit [18]. For the advective term we take the anti-symmetric spatial derivative 



d^u —7- — {u{x + a) — u{x — a)) , 
2a 



(3.2) 



where a is the lattice spacing in the spatial direction. With this choice of discretization the prob- 
lem is amenable to a local over-relaxation algorithm [p^. Starting from an initial configuration 
{u{x,t),{x,t) G A} the set of single-site variables is updated iteratively by the successive appli- 
cation of a transition probability P{u{x,t) — > u'{x,t)). We use the high-quality ranlux (pseudo) 



random number generator [20| which is essential for large-scale lattice simulations. Specific im- 



provements, e.g. Chebyshev acceleration [21] significantly reduce thermalization and autocorrela- 
tion times for the relevant observables. 

In our simulations we use periodic boundary conditions in space and fixed (Dirichlet) boundary 
conditions in time. That way we eliminate the zero mode from the dynamics. One important point 



is that the probability distribution functional ( [2.3| ) defines a stationary process for a system of 
infinite extent (in the time direction), i.e. 



{u{X[,t)u{x2,t) ■ ■ - u{x„,t)) = {u{X[,t + t')u{x2,t + t') ■ ■ ■ u{x,2,t + t')) . 



(3.3) 



In practice, this condition has to be checked explicitly. We find that for a finite space-time lattice 
this property holds to good approximation in the middle of the configurations where boundary 
effects are neglible. This defines the physical region where one may extract correlation functions. 



Another issue is that of Galilean invariance [ ]22| , |23| ]. Both the action ( |2.3| ) and the measure are 
invariant in the continuum under Galilean transformations 



x^x + r, u{x) ^ u{x + r) +v , r = vt. 



(3.4) 



To avoid an overcounting of field configurations one has to perform a gauge fixing, where one 
inserts 

(3.5) 



1 = ^[u] J dv5[u{r = vt,t)+v] 



in the dynamic functional (|2.3|). Here, ^ [u] is the Faddeev-Popov Jacobian. While gauge fixing is 
unavoidable for generic correlation functions this is not so for velocity differences that are clearly 



invariant under ( |3.4| ). 

We want to give a short remark on the computational requirements. Since we use a local 



over-relaxation algorithm, the long-range correlations imposed by the forcing ( |1.4[ ) prohibit any 
attempt to parallelize in the spatial direction. This poses a severe problem when turning to higher 
dimensions and it is absolutely necessary to switch to a global, e.g. Hybrid Monte Carlo algorithm. 
With a parallel code (in the time direction) for a 245 x 1024 lattice (space x time) our simulations 
currently run on up to 512 processors. 
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Figure 2: (a) Log-log plot of the structure function of order n ~ 5 with a linear scaling function plotted 
for comparison. Vertical bars indicate the region for the extraction of scaling exponents. Inset shows the 
local slopes versus r. (b) Structure function scaling exponents versus order n. The black curve indicates 
a bifractal scaling behavior. 



4. Results 

Structure functions are evaluated over an ensemble of configurations generated by the MCMC 
algorithm. We measure structure functions in the middle of our configurations at randomly chosen 
starting points. That way it is possible to reduce autocorrelation effects significantly. The main 
results of our simulations are shown in Fig. ^ and Fig.^. For details on the extraction of structure 
functions and the scaling spectrum we refer to |^^. Here, as an example we show a log-log plot 
of the fifth order structure function (see Fig.^). The scaling region is clearly visible, and we have 
indicated the region for the extraction of the scaling exponents by two vertical lines. In practice, 
we are bound to work at finite viscosity, and at small values of the separation we see the dissipative 
regime where the scaling breaks down. For comparison, in the inset we have plotted the local 
scaling exponents evaluated over three successive points. One may recognize, that in the scaling 
region the values lie on a plateau (as indicated by the horizontal line) which defines the scaling 
exponent. Applying this procedure to all structure functions of order n < 5 yields the scaling 
spectrum shown in Fig.^. The black line shows the bifractal scaling prediction ([L^. Our results 
are in good agreement with this prediction and also with previous results from high-resolution 



simulations employing a fast Legendre transform algorithm []13p. 

We can extract important information on the physical behavior from the probability distribu- 
tion functions (PDF) of velocity differences 3^{l^u,r). In Fig.^ we show the PDF of velocity 
differences Am = m(x + r) — u{x) for different values of the separation r plotted as a function of the 
dimensionless variable = Am/[(Am^)]'/^. One may clearly recognize the influence of the random 
forcing acting at large scales (red) where the fluctuations become Gaussian. For smaller values of 
the separation large fluctuations become strongly enhanced by the dynamics (blue, orange). In par- 
ticular, in the disspative regime (orange), at very small separations, these fluctuations are directly 
associated with the dissipative shocks (see Fig. In the intermediate range the PDF collapse (blue 
curves) and we have an indication of universal behavior. 

One particularly interesting region is indicated by the arrow in Fig. Here, the PDF for 
different values of the separation collapse exactly - this corresponds to the region |Am| <C Urms, 
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Figure 3: Probability distribution functions 3^{Au,r) as a function of the dimensionless variable = 
Am/[(Am^)]'/'^ plotted for different values of r. (a) Collapse of the PDF in the universal regime (blue). In the 
energy-containing range (red) the fluctuations become Gaussian - the random forcing dominates - whereas 
in the dissipative regime (orange) fluctuations are strongly enhanced, (b) Scaling region for the left tail of 
the PDF. The black line indicates the scahng prediction with exponent y = —4. 



r <^L where the PDF of velocity differences has the universal scaling form 

\u,r)=r-'■f(^u/r^■) 



(4.1) 



with the dynamic exponent z- In the asymptotic region —\u/r^ ^ 1 where Am < we expect the 
algebraic scaling 

^(AM,r)oc (Am)^ , (4.2) 

with 7= —4 [25, 26]. This is shown in Fig.|3|b where we have plotted the scaling region of the 
left tail of the PDFs. Though our data is not sufficient to clearly extract the scaling exponent, our 
results are in agreement with the scaling prediction. 



5. Summary 

We have demonstrated that lattice simulations can contribute to the understanding of intermit- 
tency in turbulence. Our simulations clearly show anomalous scaling for the high order moments 
of velocity differences where the exponents are in excellent agreement with previous estimates 
[1T3|]. We want to emphasize that in terms of computational efficiency our method cannot compete 
with other conventional time-advancing methods, e.g. pseudo-spectral or finite-difference methods. 
However, lattice simulations may provide a different perspective on the problem of intermittency 



where large fluctuations play a dominant role [11, 12] 
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